Evolution of Hybrid Defect Networks 



C.J.A.P. Martins 1 ' 2 '0 

1 Centro de Astrofisica, Universidade do Porto, Rua das Estrelas, 4.150-762 Porto, Portugal 
2 DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 OWA, United Kingdom 

(Dated: 14 September 2009) 

We apply a recently developed analytic model for the evolution of monopole networks to the 
case of monopoles attached to one string, usually known as hybrid networks. We discuss scaling 
solutions for both local and global hybrid networks, and also find an interesting application for 
the case of vortons. Our quantitative results agree with previous estimates in indicating that the 
hybrid networks will usually annihilate soon after the string-forming phase transition. However, we 
also show that in some specific circumstances these networks can survive considerably more than a 
Hubble time. 



I. INTRODUCTION 

The pioneering work of Kibble [l[ has shown that topo- 
logical defects necessarily form at cosmological phase 
transitions. What type of defect network is formed and 
its basic properties (such as whether or not it is long- 
lived) will depend on the characteristics of each phase 
transition, and most notably on the specific symmetry 
being broken. Understanding these processes, as well as 
the subsequent evolution of these networks, is a key as- 
pect of particle cosmology. A thorough overview of the 
subject can be found in the book by Vilenkin and Shel- 
lard @. 

Most of the work on defects in the past three decades 
has focused on cosmic strings. This is justified on the 
grounds that they are usually cosmologically benign, and 
are a generic prediction of inflationary models based on 
Grand Unified Theories [3, 4] or branes By contrast 

domain walls and monopoles are cosmologically danger- 
ous (and any models in which they arise are strongly 
constrained). However, a solid understanding of the lat- 
ter is still important, as it is the only way one can trust 
these constraints. Moreover, it is becoming increasingly 
clear, particularly in the context of models with extra di- 
mensions such as brane inflation, that hybrid defect net- 
works will often be produced. Two examples that have 
attracted considerable interest are semilocal strings and 
cosmic necklaces @, HI, HI ■ 

This is the second report on an ongoing project which 
is addressing some of these issues. In a previous work [To| 
we have developed an analytic model for the evolution 
of networks of local and global monopoles [ll|, [IH, 
The model is analogous to the velocity-dependent one- 
scale model for cosmic strings [3, llll. w hich has been 
extensively tested against field the ory |l6l. [ItJ and Goto- 
Nambu simulations [It], 13 ■ In [l(| we have also dis- 
cussed the solutions of this analytic model, and com- 
pared them with existing (relatively low-resolution) sim- 
ulations of global monopoles. Here, after a brief overview 
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of hybrid networks, we extend our analysis to the case of 
monopoles attached to one string (the so-called hybrid 
networks), and we also briefly discuss how to apply it to 
vortons. 



II. OVERVIEW OF HYBRID NETWORKS 

We start with a brief overview of previous results on 
the evolution of hybrid networks, emphasizing the dy- 
namical aspects we seek to model. A more detailed can 
be found in the Vilenkin and Shellard textbook and 
in other references that will be pointed out where appro- 
priate. We will first discuss the case of local symmetries, 
and then highlight relevant differences in the global case. 

A. Local case 

Hybrid networks of monopoles connected by strings 
will be produced, for example, in the following symmetry 
breaking scheme [3 

G^KxU(l)^K. (1) 

The first phase transition will lead to the formation of 
monopoles, and the second produces strings connecting 
monopole- antimonopole pairs; the corresponding defect 
masses will be 

Ait 

m ~ — rj m (2) 

e 

and 

M ~ Ws ■ (3) 

The characteristic monopole radius and string thickness 
are 

S m ~ (er^r 1 (4) 

6s ~ (evs)- 1 , (5) 

In the simplest models all the (Abelian) magnetic flux of 
the monopoles is confined into the strings. In this case 
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the monopoles are often dubbed 'beads'. However, in 
general (and in most realistic) models stable monopoles 
have unconfined non-Abelian magnetic charges. As in 
the case of isolated monopoles, the key difference between 
the two cases is that unconfined magnetic fluxes lead to 
Coulomb-type magnetic forces between the monopoles. 

Up to the second transition the formalism that we have 
developed for plain monopoles [Io| obviously applies, but 
the hybrid network requires a separate treatment. From 
the point of view of analytic model-building, the hybrid 
case presents one crucial difference. In the case of iso- 
lated monopoles we saw that their evolution can be di- 
vided into a pre-capture and a post-capture period: cap- 
tured monopoles effectively decouple from the rest of the 
network, and most of the radiative losses occur in the 
captured phase, where monopole-antimonopole pairs are 
bound and doomed to annihilation. This meant that the 
network's energy losses could be described as losses to 
bound pairs and we did not need to model radiative losses 
explicitly. Now compare this to the present scenario: the 
monopoles are now captured ab initio (as soon as the 
strings form) and therefore we will need to rethink the 
loss terms — as well as to account for the force the strings 
exert on the monopoles. 

The lifetime of a monopole-antimonopole pair is to a 
large extent determined by the time it takes to dissipate 
the energy stored in the string, e s ~ fiL for a string of 
length L. This is because the energy in the string is typ- 
ically larger than (or at most comparable to) the energy 
of the monopole. Monopoles are pulled by the strings 
with a force F s ~ /i ~ rj 2 , while the frictional force act- 
ing on them is F/ r , ~ 9T 2 v (where 9 is a parameter 
counting the number of degrees of freedom interacting 
with the monopoles). This friction term is of course al- 
ready included in the evolution equations of the analytic 
model [l(| EH EH : the corresponding friction lengthscale 
is if = M/(0T 2 ) ~ rj m /(6T 2 ). 

In a friction-dominated epoch the string tension should 
be compensated by the friction force, so the velocity of 
the monopoles can be estimated from 



which gives 



fi - 8T 2 v ~ 



6T 2 



(G) 



(7) 



Note that in addition to the damping due to Hubble ex- 
pansion and frictional scattering, we now have a third 
dynamical mechanism, which we may call string forcing. 
It's important to notice that in the radiation epoch the 
frictional scattering term has the same time dependence 
as the Hubble damping term, while in the matter era it 
is subdominant. 

At string formation jp^ ~ Ffrii so friction domination 
will end very quickly. This differs from the case of plain 
monopoles, whose evolution is always friction-dominated 
in the radiation era 10]. At low temperatures the string 
tension is much greater than the friction force, and a 



monopole attached to a string moves with a proper ac- 
celeration a = jjb/m ~ Tjg/r] m . 

If there are non-confined fluxes, accelerating 
monopoles can also lose energy by radiating gauge 
quanta. The rate of energy loss is expected to be given 
by the classical electromagnetism radiation formula 



()TT 



V TO/ 



JL 



(8) 



where g is the magnetic charge. This should be compared 
with the ratio of gravitational radiation losses 



-GfJ, ; 



the ratio of the two is therefore 



-grav 



'■-gauge 



mpi 



(9) 



(10) 



so if there arc unconfined fluxes the gauge radiation 
will be dominant (except in the extreme case where the 
monopoles form at the Planck scale itself). The charac- 
teristic timescales for this process to act on a monopole- 
antimonopole pair attached to a string of length L can 
therefore be written 



L 

7~rad ^ "77 



where Q 



gauge 



(ii) 



and 



(Vs/Vm) f° r g au S e radiation 
^sgrav — (Vs/mpi) 2 f° r gravitational radiation. 

These radiation losses should be included in the 
model's evolution equation for the characteristic length- 
scale L of the monopoles. It's easy to see that the corre- 
sponding term has the form 

3^ = -i^-Lf~-L = Q. (12) 

dt p 6 T ra d 

Note that in principle these terms should be velocity- 
dependent (this is discussed in [15]). However, we will 
soon see that in this case the monopoles will become 
ultra-relativistic (v ~ 1) shortly after the strings form 
and therefore the velocity-dependence is not crucial for 
the analysis. 

Another possible dynamical mechanism is that of 
string intercommuting and consequently the possibility 
of energy losses from the production of string loops. Its 
importance relative to other energy loss mechanisms is 
much harder to estimate than for ordinary strings, since 
were we're looking for an indirect effect on the monopole 
evolution, and probably one can only quantify it by us- 
ing numerical simulations. What one can say is that, as 
in the standard case, it should lead to a term in the L 
equation of the form [l5[ 



dL 
~dt 



ChybV 



(13) 



where the dimensionless coefficient Ch y b need not have 
the same value (or indeed the same order of magnitude) 
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as the standard one. We can now observe that since 
one expects to be dealing with ultra-relativistic string 
(v ~ 1) the velocity dependence will again not be crucial, 
so we can group this term with that coming from radi- 
ation losses (Eq. [P2]) and replace the coefficient Q with 
an effective which will include the effects of gauge 
radiation (if it exists), gravitational radiation and loop 
production. 

It is expected that monopole-antimonopole pairs will 
annihilate very quickly once the friction force becomes 
unimportant. The characteristic monopole velocity can 
be estimated to be approximately 



1/2 



(14) 



or v ~ I if the above is larger than unity. In the wake of 
the above discussion the corresponding rate of energy loss 
is naively estimated to be e ~ —T 2 v 2 . The approximate 
lifetime of a pair should then be 



(15) 



T 2 y 2 ' 

which in the non-relativistic case can be simply written 

m 

T nr ~ ' 

while for the ultra-relativistic case 



(16) 



J^2 



(17) 



In either case monopoles annihilate in a timescale shorter 
than a Hubble time, which in these units is 



t ■ 



mpi 

J^2 



(18) 



that there are some monopole annihilations. The strings 
become dynamically important when F s ~ F m , that is 



(21) 



at which point they pull the monopole-antimonopole 
pairs together, and the network is expected to disappear 
within a Hubble time. 

The scale-dependent monopole mass is important for 
the friction lengthscale, which is now If = M/(9T 2 ) ~ 
ij 2 n L/(9T 2 ). Another difference is that instead of gauge 
radiation we now have Goldstone boson radiation, whose 
rate of energy loss is 



tgold 



Vm 3 



(22) 



notice that this is much stronger than the energy loss 
rates due to gauge radiation (except if the monopoles and 
strings form at the same energy scale, in which case they 
will be comparable) and due to gravitational radiation 
(except if both form that the Planck scale, in which case 
they will all be comparable). 

This energy loss term should be similarly included in 
the evolution equation for the characteristic lengthscale 
L. Note that there's a crucial difference between this 
and the gauge case: here the string energy is typically 
smaller than the monopole energy, E m ~ V m L- The 
string is therefore (to some extent) irrelevant, and in this 
case the Q coefficient is Q ao id ~ cv as in the case of 
isolated global monopoles [l(|. Notice the presence of 
the velocity-dependence, although one expects that the 
velocities will always be relativistic. The effect of string 
loop production can in principle be accounted for by a 
redefined (effective) coefficient, c*. 



B. Global case 



One can also have hybrid networks of global monopoles 
connected by global strings [l9j . Just as in the case 
of plain monopoles 10], the scale-dependent monopole 
mass and the different behavior of the forces between 
monopoles in this case (long-range rather than Coulomb- 
type) changes the detailed properties of these defects and 
warrants a separate treatment. 

The tension of a global string is 

F s ~2^ s 2 ln^, (19) 

so now there's an additional logarithmic correction, while 
the long-range force between the monopoles is 

F m ~ W m ■ (20) 

If T] m >> rj s the monopoles initially evolve as if they 
were free, with L ~ t. Note that this scaling law implies 



III. MODELING HYBRIDS 

The above overview leads to the expectation that the 
network will annihilate shortly after the strings form. In 
fact we will see below that there are circumstances where 
the network can survive considerably more than a Hubble 
time, although it has to be said that we do expect this 
to be the exception rather than the rule. 

We will therefore proceed fairly quickly, discussing only 
the late-time regime where the force due to the strings is 
dominating and the hybrid networks are therefore about 
to annihilate. As discussed in Sect. |TT1 in the local 
case this will happen very shortly after the string- forming 
phase transition. In the global case, however, there will 
be an intermediate epoch where the monopoles evolve as 
if free (because the strings are comparatively very light), 
and the strings only become dynamically important at 
an epoch given by Eq. (|2"Tj) . 
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A. Local case 

In this case our evolution equations for the character- 
istic lengthscale L and RMS velocity v of the monopoles 
take the form 



dL 
3 — 

dt 



3HL + v L H 



9T 2 \ 

Vr, 



(23) 



where includes the energy loss terms from gauge radi- 
ation (if it exists), gravitational radiation and loop pro- 
duction discussed earlier, possibly with some coefficient 
of order unity, and 



dv 
~dl 



L 



V, 



„L 2 d 



H 



H 



OT 2 



(24) 

The velocity equation now has two accelerating terms, 
due to the strings and the inter-monopole Coulomb 
forces, which we parametrize with coefficients k s and k m 
which are expected to be of order unity. An exception is 
the Abelian case, where there are no Coulomb forces, so 
k m = in this case. 

It's important to realize that friction will play a crucial 
role in the radiation era, where it is more important than 
the Hubble damping term (the opposite is true for the 
matter era). Indeed in the radiation era we can write 



f 



T 2 



= H + e — = 



i / 1 



} m P i 

]rn 



hi 
t 



(25) 



where in the last step we have defined an effective coef- 
ficient that is usually much larger than unity (except if 
there were no particles interacting with the string, 8 = 0). 
On the other hand, in other cosmological epochs (in par- 
ticular the matter-dominated era) the friction term is 
negligible, and we have A* = A (where we are defining 
a oc t x ) as usual. 

It is illuminating to start by comparing the various 
terms in the velocity equation. The ratio of the monopolc 
and string accelerating terms is 



F, 



1 



k s LdHVs 



(26) 



which is always much less than unity (except if A; s hap- 
pened to be extremely small). The Coulomb forces are 
always negligible relative to the string forces. This was 
expected, since as we pointed out earlier the energy in the 
strings is typically larger than that in the monopoles. A 
useful consequence is that we do not need to treat the 
Abelian and non- Abelian cases separately. 

Now let us compare the damping and string accelera- 
tion terms 



Fd 



A* 



Vn 



mpi 



(27) 



where in the last step we assumed that we are in the 
radiation era. Given that the evolution of the monopoles 



before string formation leads to a scaling law v oc a -1 for 
monopole velocities [10(, which are therefore extremely 
small when the strings form, the above ratio is always 
less than unity. 

This analysis therefore quantitatively confirms the 
naive expectation that as soon as the strings are formed 
the string acceleration term will dominate the dynamics 
and drive the monopole velocity to unity. Recalling that 
the initial monopole velocity can effectively be taken to 
be zero, we can write the following approximate solution 
of the monopole velocity equation 



In- 



2f s (t-t s ) 



(28) 



and we can compare the epoch at which the monopoles 
become relativistic (t c ) with that of the string-forming 
phase-transition (t s ), finding 



t s mpi 



(29) 



so they become relativistic less than a Hubble time af- 
ter the epoch of string formation; notice that this ratio 
depends only on the energy scale of the monopoles, and 
not on that of strings. 

We can now proceed to look for solutions for the char- 
acteristic monopole lengthscale L, assuming for simplic- 
ity that v = 1 throughout. We can look for solutions of 
the form L oc t a for generic expansion rates (a oc t x ), and 
it is straightforward to see that there are two possibilities. 
There is a standard linear scaling solution 



Q* 



3(1 - A) - A, 



t 



(30) 



in which the energy loss terms effectively dominate the 
dynamics. In this case the monopole density decays 
slowly relative to the background density, 



Pm ,_i 
— oct 

Pb 



(31) 



This scaling law can in principle exist for any cosmolog- 
ical epoch provided A < 3/4 (for example, in the matter 
era we have L = 3Q*t). However, in the radiation epoch 
we would need the unrealistic 9mpijr\ m < 2, which effec- 
tively would mean that friction is absent (9 = 0). 

The alternative scaling solution, for epochs when fric- 
tion dominates over Hubble damping (such as the ra- 
diation era) and also for any epoch with A > 3/4 (even 
without friction) has L growing with a power a > 1 given 
by 



a = A + - A* = - I 4A 
3 3 1 



(32) 



if there is no friction the scaling law can simply be written 
L oc a 4/3 , (33) 
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from which one sees that although the Hubble term is 
dominant, the scaling is faster than conformal stretching 
(L oc a) because the velocities are ultra-relativistic. The 
important point here is that in this regime L grows faster 
than L oc t, and the number of monopoles per Hubble 
volume correspondingly decreases. 

At the phenomenological level of our one-scale model, 
this corresponds to the annihilation and disappearance 
of the monopole network. We can easily estimate the 
timescale for this annihilation — it will occur when the 
number of monopoles per Hubble volume drops below 
unity (or equivalently L > dn)- Assuming a lengthscale 
L s = st s at the epoch of string formation (note that the 
evolution of the monopole network before the stings form 
is such that s can be much smaller than unity) , we easily 
find 



ts 



1 



3 (2 
9 \s 



1 



Vn 



mpi 



(34) 



which is comparable to the estimate we obtained using 
the velocity equation, Eq. (f2"5)) . Notice that this is a 
much faster timescale than the one associated with radia- 
tive losses, which can therefore be consistently neglected 
in this case. There's nothing unphysical about this 'su- 
perluminaP behaviour, as explained in (20j . We do have a 
physical constraint that the timescale for the monopole 
disappearance should not be smaller than the (initial) 
length of the string segments, t a > L s , but this should 
always be the case for the above solutions. 

This analysis shows that the monopoles must annihi- 
late during the radiation epoch, if one wants to solve the 
monopole problem by invoking nothing but a subsequent 
string-forming phase transition. Any monopoles that sur- 
vived into the matter era would probably be around to- 
day. It also shows that monopoles can also be diluted 
by a sufficiently fast expansion period. Inflation is of 
course a trivial example of this, but even a slower expan- 
sion rate 3/4 < A < 1 would be sufficient, provided it is 
long enough to prevent the monopoles from coming back 
inside the horizon by the present time. 



The evolution of a given segment of a local string of 
length I in the context of the VOS model has been stud- 
ied in 21]. In the present notation the evolution equation 
has the form 



df 
~dt 



(1 - 2v 2 )H£- 



-v 2 - 



(35) 



note that as expected energy loss mechanisms reduce the 
length of the string segment. There's an analogous equa- 
tion for the string segment's velocity, but we already 
know we can safely assume v ~ 1. Using this and the 
specific form of the friction term, the above equation can 
in fact be written 



d£ = _A* 
dt t 



(36) 



which can easily be integrated to yield (for I = £ s = st s 
at the string formation epoch t = t s ) 



l + A* 



(37) 



We emphasize that for simplicity we are assuming that 
all string segments are formed with the same length: at 
a detailed level this is unrealistic, but it is nevertheless 
sufficient to provide reliable qualitative estimates. 

We can now estimate the monopole annihilation epoch 
by simply looking for £ ~ 0. As before, the answer will 
depend on the relative importance of the the radiative 
loss and friction terms. If friction is negligible and the 
radiative losses dominate, then we find the radiative de- 
cay timescale 



tr f. , 1 + AA V1+A * 



(38) 



which is much larger than unity. On the other hand, if 
friction is important and the Q term is subdominant the 
annihilation timescale is the much faster 



B. An aside: string evolution 

Notice that in the above we haven't yet said anything 
about the evolution of the strings. One could try to de- 
fine a string characteristic lengthscale in the usual way, 
Ps = AW ' L 2 S . However, one should be careful about doing 
this, since in this case there is no a priori expectation 
that the distribution of string segments connecting the 
pairs of monopoles and antimonopoles will form an ef- 
fectively Brownian network. (In fact it is quite likely 
that it doesn't, although this is something that warrants 
numerical investigation.) 

A safer and simpler way of describing them is to look 
at the evolution of an individual string segment, and 
then make use of the fact that we already know how the 
monopole number and characteristic lengthscale evolve. 



h. „ i + -3i 

t.*t 



mpi 



(39) 



which is again comparable to our previous estimates, 
given by Eqs. and ([M]) . 

Since each piece of string is connecting two monopoles, 
a very simple estimate of the total length in string in a 
Hubble volume (and hence of the string density) is ob- 
tained multiplying half the number of monopoles in that 
volume by the typical length of each segment. This leads 
to 



Ps 



Ik 



At) 



(40) 



and as expected the string density decays relative to that 
of the monopoles. 
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Global case 



In the global case the evolution equations will be [l(| 



dL 



L 



3—=3HL + v z - + c*v 

at id 



(41) 



and 

dv 
~dt 



(l-« 2 ) 



L \ 3/2 + « ln ^ 
L ^ 2 „ 5 S 



and the scale dependence of the monopole mass and 
string tension imply that this case differs in two ways 
from the local case. 

Firstly, the damping term at the string formation 
epoch (assumed to be in the radiation epoch) can be 
written 



(43) 



The friction term is now subdominant relative to Hub- 
ble damping, and of course it decreases faster than it. 
Friction can therefore be ignored in the analysis. 

Secondly, the monopole acceleration due to the 
(global) strings now has an extra logarithmic correc- 
tion, but more importantly the logarithmically divergent 
monopole mass again implies that this force is inversely 
proportional to L (as opposed to being constant in the 
local case). The large lengthscale inside the logarithm 
should be the string lengthscale, but we have substituted 
it for the monopole one, which we'll simple denote L. 

Starting again with the velocity equation, the ratio of 
the string and monopole acceleration terms is now 



t-1 t 2 rfi/ 1 ^ 



■ In (Lr) s ) 



(44) 



and in particular at the epoch of string formation we have 



- 1 - in 



km Vrn Vs 



(45) 



so now F s is initially sub-dominant, except if we happen 
to have rj s ~ rj m . Moreover, given that L « t with a 
proportionality factor not much smaller than unity, while 
monopoles evolve freely (before the effect of the strings is 
important), the ratio will only grow logarithmically, and 
so the effect of the strings might not be felt for a very 
long time. Specifically this should happen at an epoch 



Vs I Vm 

mpi \ <qj 



(46) 



naturally this is exactly the same as Eqn. (f2lj) . 

While strings are dynamically unimportant we have 
v = const, as for free global monopoles, and even when 
they become important the velocity will still grow very 



slowly (logarithmically) towards unity. Hence, although 
the ultimate asymptotic result is the same in both cases 
(v — > 1), the timescale involved should be much larger in 
the global case. 

Moreover, recall from [lCj that in the gauge case the 
monopole velocities before string formation were neces- 
sarily non-relativistic and indeed extremely small, and it 
is the strings that make them reach relativistic speeds. 
In the global case this is usually not so, as the monopoles 
will typically have significant velocities while they are free 
(although their magnitude depends on model parameters 
that need to be determined numerically). We can there- 
fore say that the impact of the forces due to the strings 
is vastly smaller in the global case. All this is due to the 
fact that the force due to the strings is inversely propor- 
tional to L instead of being constant. 

As for the evolution of L, at early times (before strings 
become important) we have L oc t just like for free global 
monopoles. Eventually the strings push the monopole 
velocities close(r) to unity, and in this limit the L evolu- 
tion equation looks just like that for the global case with 
the particular choices A* = A and Q± — c. However, we 
must be careful about timescales, since here the approach 
to v = 1 is only asymptotic. 

Bearing this in mind, for A < 3/4 we will still have 
linear scaling solution 



L = 



3-4A 



(47) 



notice that this is exactly the ultra-relativistic (v = 1) 
scaling solution we have described in [l(| — cf. Eqs. (70- 
71) therein. Therefore, if the network happened to be 
evolving in the other (subluminal) linear scaling solu- 
tion, the only role of the strings would be to gradually 
switch the evolution to the ultra-relativistic branch. This 
scaling law can in principle occur both in the radiation 
and in the matter eras, and it follows that in this case 
the monopoles will not disappear at all, but will continue 
to scale indefinitely (with a constant number per Hubble 
volume). 

In this case an analysis in terms of the length of each 
string segment would have to take into account an initial 
distribution of lengths. Moreover, since the initial ve- 
locities need not be ultrarelativistc, the larger segments 
(which will have smaller coherent velocities) should grow 
at early times, while the smaller ones will shrink. The 
decay time will obviously depend on the initial size. We 
note that such a behavior has been seen in numerical 
simulations of semilocal strings [22| . 

The alternative scaling solution, for any epoch with 
A > 3/4 has L growing with as 



L oc a 4 / 3 , 



(48) 



and again corresponds to the annihilation and disappear- 
ance of the monopole network, which would occur at 



t„ 



1 



tr [(1 - A)rf/( 4A - 3 > 



(49) 
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for an initial lengthscale L s — rt r ; given that we now ex- 
pect r to be not much smaller than unity, this is likely to 
be very soon after the onset of this scaling regime — that 
is, very soon after the velocities become ultra-relativistic. 
(Notice that in the local case the monopoles became 
ultra-relativistic very soon after the strings formed so 
t c ~ t s , but in the present case t c » t s , except if the free 
monopoles were already evolving in the ultra-relativistic 
branch.) 



IV. TOPOLOGICAL STABILITY 

There is in principle a further energy loss mechanism 
that we have not discussed so far. These strings are 
not topologically stable: they can break, producing a 
monopole-antimonopole pair at the new ends. This pro- 
cess was first discussed in [l9[ and more recently in [23j . 

This is a tunneling process (usually called the 
Schwinger process), and its probability per unit string 
leng th and per unit time has been estimated to be 



a I m 

P ~ exp — 7T 

2?r V M 



(50) 



By substituting Eqs.((2HHJ) onc notices that in typi- 
cal GUT models the exponent is very large, say ~ 
10 3 (?7 m /»7a) 2 . This implies that this probability is neg- 
ligibly small, even in the most favorable case where 
the strings and monopoles form at the same energy 
scale. Therefore, for all practical purposes such strings 
(whether they are local or global) can be considered sta- 
ble. 

Nevertheless, it is conceivable that in string-inspired 
models such as brane inflation there are models (or re- 
gions of parameter space in some models) where the ex- 
ponent is not much larger than unity. Here we briefly 
discuss what happens in this case. A different scenario, 
where the bead- and string-forming phase transitions are 
separated by an inflationary epoch (and thus the beads 
are initially outside the horizon), has been recently dis- 
cussed in [23j | . 

Again it is simplest to look at the effect of this addi- 
tional energy loss term on an individual string segment. 
It's straightforward to show that Eq. (f3"6")> now becomes 



(51) 



at t 



Assuming that the exponent in the Schwinger term is 
of order unity, it's easy to see that for large segments 
(that is those much larger than the string thickness S s ) 
the Schwinger term is the dominant one, and it quickly 
breaks up the segment into smaller ones. On average 
(and neglecting the standard damping and energy loss 
terms) the length of a given segment is expected to shrink 



and the epoch (relative to that of string formation) at 
which this process would, on average, have reduced the 
segment to a size 8 S is 



t s e mpi 



(53) 



which is typically less than a Hubble time. As the seg- 
ments shrink the standard damping and especially the 
radiative energy loss terms will become more important, 
and we would eventually switch to a solution of the type 
discussed earlier, cf. Eq. (|37[) . The segment size when 
this happens will be very small, so s « 1 in Eq. (j38j) . 
and these segments should therefore decay in much less 
than a Hubble time. 



V. A DIGRESSION: VORTONS 

Since vortons [26| are to a large extent point-like ob- 
jects decoupled from the string network that produced 
them, it would be interesting to use the monopole model 
to describe the evolution of their overall density. 

The most naive approach would be simply to say that 
there are no inter-vorton or other forces (other than fric- 
tion and Hubble damping) affecting their dynamics. This 
immediately implies that there will be no acceleration 
term in the velocity equation (k = 0), and hence vortons 
will necessarily be non-relativistic. 

At the same naive level there should be no energy 
losses to account for in the lengthscale evolution equa- 
tion. Classically this is justifiable by saying that we 
can assume that all vortons form around the same epoch 
(related to the superconducting phase transition, which 
need not be the same as the string-forming phase transi- 
tion) and that once formed a vorton stays there forever 
(or at least has a very long lifetime) — these issues are 
discussed in [2|. Of course we're not addressing the is- 
sue of the vorton density at formation (which would set 
the initial conditions for the model). A different issue 
which we are also neglecting is that there could well be 
an energy loss term due to quantum effects. 

With these caveats in mind, it follows from our previ- 
ous work [T3| that the evolution equations will be 



and 



3 f = (3 + -^ + ^ 



(54) 



(55) 



Note that here L is a characteristic vorton lengthscale, 
which is related to the vorton density through 



Pv 



L 3 ~ L 3 



(56) 



£(t) 



l + P£.(t-t„) 



(52) 



where rj s is the string symmetry breaking scale (/x ~ rjl) 
and £ v is an effective vorton length (not necessarily its 
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radius/size, since there will be a significant amount of 
energy associated with the charge/current). For simplic- 
ity (and consistently with the above approximations) we 
will further assume that this length is a constant. 

The simplest possibility we can consider is the case 
where friction is negligible (if — > oo). The analysis of 
this case is trivial, and it leads to 



T 1 1 

L cx a , v (X — oc — 
a L 



(57) 



This of course leads to p v cx a~ 3 , which is in fact what has 
been found in previous analysis both for generic vortons 
[27l ] and for the specific case of chiral vortons [Hj]. In 
either case there are several possibilities for the initial 
density, depending on details of the underlying particle 
physics model. There is no published analysis of vorton 
scaling velocities, although the effect of superconducting 
currents on the loop velocities has been discussed in [2^] , 
and the above scaling law for velocities is well known for 
defects in condensed matter systems. 

However it's also important to consider the role of fric- 
tion, which in this context should not have the same form 
as the standard case. This is discussed in [3(| and also 
in the Vilenkin and Shellard textbook 0] • Because a su- 
perconducting string creates a magnetic field 



B s (r) ~ 2- 
r 



(58) 



(J being the current) a string moving in a plasma is 
shielded by a magnetocylinder which contains the mag- 
netic field which does not allow the plasma to penetrate. 
The magnetocylinder radius is 



J 



v VPp 



(59) 



where p p is the plasma density, and this leads to a friction 
force per unit length 



F p ~ JVyfpp~ 



(60) 



which in our notation corresponds to a friction length- 
scale 



Note that the current can have a maximum value 



■J I, 



(61) 



(62) 



We can now study the effect of this friction term on 
the above evolution equations. It's easy to see that the 
scaling of the characteristic vorton scale is unchanged, in 
other words we still have L cx a. However, there is an 
effect on the velocity equation, which now takes the form 



dv 2 



H 



(63) 



Now for a universe with a critical density 

3H 2 



Pp 



8ttG 



(64) 



so in fact the two terms are proportional to the Hubble 
parameter, and which of the two dominates will depend 
on the strength of the current J. For a generic expansion 
law a cx t x we can write the scaling law for the vorton 
velocities as 



v oc t' 



where 



= X 1 + 



-0 



3 Jmpi 



87T 7?2 



(65) 



(66) 



For small currents we have v cx t~ x cx a -1 as before, but 
the effect of the plasma is to reduce the velocities at a 
much faster rate. Note that for a maximal current the 
scaling exponent will be approximately (mpi/r) s ), which 
is much larger than unity. 

It would be interesting to test this behavior numeri- 
cally, as well as to understand the magnitude and form 
of energy loss terms due to quantum effects. 



VI. SUMMARY AND OUTLOOK 

We have applied a recently developed analytic model 
for the evolution of monopole networks [l(| to the case 
hybrid networks of monopoles attached to a single string, 
and also to the case of vortons. We discussed scaling so- 
lutions for both local and global hybrid networks, gener- 
ically confirming the expectation that the network will 
annihilate shortly after the strings form, but also high- 
lighting the fact that there are circumstances where the 
network can be long-lived — although it has to be said 
that we do expect this to be the exception rather than 
the rule. 

Studying the evolution of hybrid defects may seem 
relatively uninteresting because of their brief existence 
or even due to the topological instability discussed in 
Sect. IIV1 However, this study has a second moti- 
vation. We have studied local strings attached to lo- 
cal monopoles and global strings attached to global 
monopoles. These are two simple test cases before tack- 
ling a third, more interesting case: local strings attached 
to global monopoles. This corresponds (modulo a few 
subtleties) to the case of semilocal strings [3l|, |32j , which 
we will address in a subsequent publication. 
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